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Abstract. We investigate the exact dynamics of a system of two independent 
harmonic oscillators coupled through their angular momentum. The exact analytic 
solution of the equations of motion for the field operators is derived, and the conditions 
for dynamical stability are obtained. As application, we examine the emergence of 
squeezing and mode entanglement for an arbitrary separable coherent initial state. 
It is shown that close to instability, the system develops considerable entanglement, 
which is accompanied with simultaneous squeezing in the coordinate of one oscillator 
and the momentum of the other oscillator. In contrast, for weak coupling away from 
instability, the generated entanglement is small, with weak alternating squeezing in 
the coordinate and momentum of each oscillator. Approximate expressions describing 
these regimes are also provided. 


Exact dynamics and squeezing in two harmonic modes coupled through angular momentum2 

1. Introduction 

Models based on coupled harmonic oscillators have long attracted attention in several 
different fields due to their wide range of applications [H [2] El IH [5j El [H 8j|9] [lOj. In 
particular, the case of two harmonic modes coupled through their angular momentum, 
which describes the motion of a charged particle within a general harmonic trap in a 
uniform magnetic field or, equivalently, the motion in a rotating anisotropic harmonic 
potential [HI H2 H3], has been employed in distinct scenarios, such as rotating nuclei 
[13], quantum dots in a magnetic held p4j and fast rotating Bose-Einstein condensates 
mm m gej within the lowest Landau level approximation [El can EH]. Since its 
Hamiltonian is quadratic in the held operators, the model is also suitable for simulation 
with optical techniques [22] . 

In a previous work [23] its dynamics in coordinate representation was analyzed in 
detail, showing that it exhibits a complex dynamical phase diagram, with stable as 
well as distinct types of unstable (i.e., unbounded) dynamics. We have also examined 
the generated entanglement between the modes, both in ground and thermal states 
[24 ] (vacuum and thermal entanglement) as well as that obtained after starting from 
a separable vacuum state [25]. It was shown, in particular, that this system is able to 
mimic typical entanglement growth regimes arising after a quantum quench in complex 
many body scenarios [26]. Entanglement is of course essential for quantum information 
applications [27], and a large entanglement growth with time after starting from a 
separable state in a many-body system, is indicative of a system dynamics which cannot 
be efficiently simulated by classical means. 

On the other hand, quantum squeezing constitutes another topic of great current 
interest PSIESIEDIEI]. Its relation with entanglement has been investigated in different 
systems [28l 29] [30] 31) [32] 33] [34], with entanglement normally inducing squeezing in 
certain observables. In particular, in [32] the exact dynamics of entanglement and 
squeezing in a two-mode Bose-Einstein condensate interacting through a Josephson- 
like coupling was determined. Squeezing is important for quantum metrology, i.e., for 
improving the accuracy in quantum measurements [3T1 [35] . and it has been shown that 
in some cases spin squeezing can be employed to detect entanglement [28) [30, [31] [33]. 
Nonetheless, in some regimes (like the linear case of the model considered in [32]) 
squeezing may also arise without substantial entanglement. 

In this work we first derive the exact analytic expressions for the temporal evolution 
of the Heisenberg held operators of two harmonic modes coupled through their angular 
momentum. The obtained result is valid for all values of the system parameters, i.e., in 
stable as well as unstable dynamical regimes, and allows to determine the exact evolution 
of an arbitrary observable of the system. We then apply this result to determine 
and examine the dynamics of squeezing, which has so far not been investigated in 
this model, and its relation with the generated entanglement, when starting from a 
separable coherent initial state. We will show that different regimes can arise depending 
on the value of the rotational frequency. Close to the instability point, appreciable 
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entanglement is generated, accompanied with simultaneous squeezing in one of the 
variables of each mode, while for small couplings, the generated entanglement is weak, 
with small squeezing appearing in both variables of the mode at alternating times. 
Approximate simple expressions describing these two distinct regimes are also provided. 


2. Formalism 


2.1. The model 


The Hamiltonian of the system under study can be written as 




rmjj\ 1 n 2 muj l ^ 

- -Qi H- P-2 H- -Q 

2 2 m 2 2 v 


UJ {Q\P 2 — P1Q2) , 


( 1 ) 


where the subscripts 1 and 2 refer obviously to the first and second oscillator. The 
oscillator frequencies uq, cu 2 and the rotation frequency 00 will be assumed real and 
satisfying, without loss of generality, 

oji > co 2 > 0 , 00 > 0 . (2) 


In the case of a particle of charge e in a magnetic field H along the z axis within a 
harmonic trap with spring constants Ad, K2 in the x, y plane, we have co = e\H\/( 2 mc), 
with uij = Kj/m + co 2 (j = 1,2) [23] (motion along z is obviously decoupled from 
that in x,y plane, which is that described by (fill ). Eq. (OQ) also represents the intrinsic 
Hamiltonian describing the motion (in the x, y plane) of a particle in a harmonic trap 
with constants mco 2 rotating with frequency c0 around the z axis [23] . 

By expressing the position (Qj) and momentum (Pj) operators in terms of the 
dimensionless annihilation (aj) and creation (aj) boson operators, Qj = + a j)j 

Pj = —iJ , PP^{aj — Oj), j — 1, 2, we can rewrite (H]) as 


where 


H — hoj\ ^ + hoj-2 + 7 ^ 

— ih\- (aia .2 — , 



A ± — oj 


f Ui ± u 2 \ 
\2y/U\U2 J 


( 3 ) 

( 4 ) 

( 5 ) 


In the isotropic case oj\ = ui 2 = cuo, A_ = 0 and both H and the angular momentum, 
which becomes just —ih\ + {a\ai — a|a 2 ), commute with the total boson number N = 
ajai + a\a 2 and between themselves. 

However, in the anisotropic case 7 - 002 , the angular momentum term does not 
commute with H and no longer conserves the total boson number. This entails, in 
particular, that in contrast with the isotropic case, the system vacuum will become 
entangled as co increases, and that an initially separable state | 0 i)| 0 2 ) (product of the 
vacuum of each oscillator) will become entangled as time increases if co 7 ^ 0. In addition, 
the system may become unstable if co increases sufficiently, as discussed below [23]. Note 
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that the exact dynamics of the special case oj\ = u> 2 with a Q 1 Q 2 coupling, which was 
studied in detail in [ 4 ], corresponds formally to the case A + = A_, after replacing 
a 1 —> ia±. Also, the two-mode Bose Einstein condensate model of [32] (which contains 
as well non-linear terms) would formally correspond in the linear case to A_ = 0 and a 
time-dependent A+. 


2.2. Exact solution 


Let us now derive the explicit solution of the Heisenberg equations of motion for the 
held operators, 

ihdj = [aj,H] , J = 1,2. (6) 

Eq. ([6]) leads to the linear system 

d\ — —ictqcq T A_|_(Z2 T A_o)> 
d 2 = — ico 2 a 2 — A + ai -|- A_oj 
which can be written in matrix form as 


( 7 ) 


a 

at 


= n 


a 


a 1 


( 8 ) 


where a = 


a 1 
a 2 


at = 


j and E is the 4x4 non-hermitian matrix 


U = 


( w>i i\ + 0 iA_ \ 

— i\ + co 2 iX- 0 

0 iA_ — oji iX + 

y iX- 0 — iX + —u> 2 J 

The exact solution of Eq. (JBJ) can be expressed as 


( 9 ) 


o(() 

o*(«) 


=m 


0 ( 0 ) 

a+(0) 


where 


U(t) = exp [—iHt\ = 


U(t) V(t) 

V‘(t) U‘(t) 

is a 4 x 4 matrix satisfying (I denotes the 2x2 identity matrix) 


U(t)MU\t) =M, M = 


I 0 
0 -I 


( 10 ) 


( 11 ) 


( 12 ) 


since = E. Eq. ensures the preservation of the equal time commutation 

relations V t\ 

[Oi(t ), a](t)} = [UU^-VV%(t) = S tj , K(t), aj (t)] = [UV'-VU*]^) = 0 , (13) 

implying that Eq. (fTOl) represents a proper time-dependent Bogoliubov transformation 
for the held operators. 
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Setting in what follows a,j( 0) = aj, at(0) = at, Eq. (HU leads explicitly to 

a j(t) — Uji(t ) a\ + Uj 2 (t ) a-2 + Vji{t) o\ + Vj 2 {t) a\ , ( 14 ) 


and the corresponding adjoint equations for at(£), where the elements Ujk(t), Vjk(t ) can 
be obtained from Eq. (ITT]) through the diagonalization of EL: 


with 




Vjj{t) 


u 12 (t) 

v 12 (t) 


U 21 (t) 


= -{(1 4" 7a) COSUJ +t + (1 — 7 j) cosca_f 

. r ._, ...sinca+f _ _ . sincaA... 

- iuj[{ 1 + 8j )——-f (1 — 8j )———]} , 




= i 


^ 7 - +1 caia; 2 A + A sincu+t 


ujjA 


LU- 

sinca_t, 

C0- 1 


(15) 

(16) 


_ A± (cai A oa 2 ) 2 sin ta + f (ag ± ta 2 ) 2 sintaA 

2 ^ 2A ca + [ 2A ca_ 

. .Vl±U2, . ,xi 

+ i -x-(cosca+r — cosca_r)J, 


A 

-U 12 (t), V 21 (t) = V? 2 (t), 


UJn 


7j_l J 2A ’ J 7j+ 2Aw? 


(17) 

(18) 

(19) 


Here ca± are the system eigen frequencies, i.e., the eigenvalues of the matrix EL (which 
are ±ca + , ±ca_), given by 


~ A 2 + ±A-y/c a 2 + ^ ± A , 

(20) 

where 


A — ^/a + (u;i + ca 2 ) + (cai ca 2 ) [ 4 A_] 

(21) 

- + 2W 2 (W 2 + W 2) . 

(22) 

It is verified that for ca = 0 (A± = 0), A = AAA and hence ca+(_) 

= cai( 2 ), leading to 


Ujk(t ) = 8jke~ Wjt , Vjk(t) = 0. The free evolution cy (t) = e~ lUjt aj is then recovered. 
Also, in the special case cai = ca 2 = ca and A + = A_ = k we recover the eigenfrequencies 
U!± = \/uj 2 ± 2/ttU of H]. 

On the other hand, in the isotropic case u\ = ca 2 = we have A_ = 0, implying 
A = 2caca 0 and 


w± = ca 0 ± ca , H(t) = 0 , (cai = ca 2 = ca 0 ), (23) 

which leads finally to 

ai(t) = e~ luot {ai coscaf+a 2 sin cat), a 2 (t) = e _ *‘ J ° t (—ai sin cat+a 2 coscat) .(24) 

This is equivalent to a beam-splitter type transformation of angle ut of the field 
operators. In this case the angular momentum term commutes with H and just rotates 
the field operators with angular frequency ca. 
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The general solution (I14p - (|18p can be derived by many other methods. For instance, 
we may write the solution of Eq. TO for each operator as 

a j(t) = e lHt/h a je - tHt/h = a, + j [H, a,] + ^ (^j [H, [H, aj ]] + ... (25) 

which leads immediately to the form (fl4P for a,j(t). And insertion of a trial solution of 
the form (1141) in (J?P (i.e., the so called Sen-Mandal approach [MUST]) leads to a linear 
system of first order differential equations for the coefficients Uj k (t), Vj k (t), namely 
iU = WJ , i.e., 


U k 

v; 


= n 


u k 

v; 


, A: = 1,2 


(26) 


where Ui, = 


are the k th column of U, V. Eq. (051) and 


U lk \ / V lk 

u 2k y k \ v 2k 

the initial conditions Uj k (0) = Sj k , Vj k ( 0) = 0 lead again to the solution (115]) — (TTSP . We 
finally notice that the system (J7|) leads, after successive derivation, to the fully decoupled 
quartic equations 

aj +(cj^+cj 2+2(A^—A 2 ,))dj+((A_|_+A_)^— cjicJ2)((A+—A_) 2 — u)\UJ2)aj = 0(27) 

for j = 1,2, which lead at once to the eigenfrequencies (1201) (after inserting a trial 
solution aj(t) oc e iat ) and again to the solution (ITsp - flTsp after inserting the initial 
conditions for the operators and their derivatives. The matrix U(t), and hence all 
coefficients Uj k (t ), Vj k (t), also satisfy Eq. (l27p . 


2.3. Dynamical stability and normal mode decomposition 

In the general case, a close inspection of the eigenvalues (1201) reveals that u j± are both 
real and non-zero only if A + + A_ < or A + — A_ > i.e., if a; < a; 2 or 

u> > uii, which is equivalent to 

(uj — — u} 2 ) > 0 ■ (28) 

Dynamical stability (bounded quasiperiodic dynamics) is then ensured if Eq. Q2HP is 
satisfied. On the other hand, if A + — A_ < yJuj\UJ 2 < A+ + A_, i.e., co 2 < uj < uji or in 
general 

{uj — — oj 2 ) < 0 , (29) 

oj. |_ remains real but io_ becomes imaginary (cu_ = i|o;_|) implying that the dynamics 
becomes unbounded. In this case we should just replace 
sino;_f 


-4 


0J- 


sinli 

a; 

— 

t 


\u 




COS UJ_t -4 cosll \uj-\t 


(30) 


in Eqs. (TTSP — (jT8l) . entailing that all operators “increase” (i.e., deviate from their initial 
values) exponentially with time. Nevertheless, Eq. (TT2l) and hence the commutation 
relations (TT3P remain satisfied. 
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Finally, if co = co 2 or co = (Ui, i.e., if 


(u - ui)(u - u 2 ) = 0, (31) 

we obtain a critical regime where co + > 0 but co- = 0, in which case the matrix E 
becomes non-diagonalizable unless co i = co 2 (Landau case). If co_ = 0, we should just 
replace the corresponding expressions by their natural limits, i.e., 
sinu;_£ _ 

- —y t ? cos (jj—t — y 1, (32) 


in Eqs. (TT5|) — (jT8|) (the ensuing solution also follows from (HI]) after using the Jordan 
canonical form of E for cu_ = 0 |23]), entailing that the dynamics is again unbounded 
if co\ ^ co 2) with the deviation from the initial values increasing now linearly with time. 
Eqs. (TT21) are again preserved. In the Landau case co = coi = co 2l the coefficients of the 
linearly increasing terms vanish and the dynamics is again bounded, given by Eq. (123|) 
for co = co 0 - 

A standard normal mode decomposition of H becomes feasible in the dynamically 
stable phases (co± > 0) [23]. In terms of the standard normal mode boson operators 
6 ±, given in the appendix, we may express (EQ) in the first dynamically stable sector 
co < co 2 as 


H = kjO + (b\b + + \) + Tuo_(tf_b_ + |), co < co 2 (33) 

with b±(t ) = e~ lul±t b±(0 ), whereas in the second dynamically stable sector co > uo\ we 
have 


H = hio + (b\b + + |) — hcn_(6t&_ + |), co > coi , (34) 

with b_(t) = e lu ~ t b^( 0) (and o;_ > 0). Eq. (IH4|) entails that in this region, the system is 
no longer energetically stable. 


3. Application 

3.1. Squeezing and entanglement 

We have now all the elements for investigating the evolution of distinct quantum 
properties of the system, such as entanglement and squeezing. We start by noting 
that the number operators for each mode are given by (here j , k, l — 1, 2) 

Nj(t) = a]{t)aj(t) = [Uj k (t) Uji (t) a\ai + E/ fc (t)E j7 (t)a fc af 

k,l 

+ U* k (t)Vji(t)alal + V* k (t)Uji(t)a k ai \, (35) 

indicating that they will acquire a non-zero average even if there are initially no bosons: 
If the system starts at the separable vacuum 100) = |0i)|0 2 ), where OjjOj) = 0, from 
Eqs. (1551) and (fTH - (fT5jl we obtain, setting ( O) 0 = (0010100) and j, k = 1, 2, 

(W))»= Ei^wi 2 

k 


( 36 ) 
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0J 2 {pJi - OJ 2 )' 


£ 

v=± 


+ 


161 WIUJ 2 
U) 2 (UJ1 + LU 2 ) 2 UJiUJ2 


(1 


. (OJI — U2) 2 \ sintivf 
2A / 


+ iv 


(uJl — U)2 ) cos UJ u t 


sin lj -(-1 

UJ + 


sine j—t 


(37) 


which will be normally non-zero for t > 0 unless u)\ — u; 2 . ( N 3 (t)) is then proportional 
to (cui — cu 2 ) 2 , and is larger in the mode with the lowest frequency ojj (due to the last 
term in (1371)1. 

If the initial state is instead a product |aqa 2 ) = \oli)\ql 2 ) of coherent states | a 3 ) for 
each oscillator, with a 3 \a 3 ) = a 3 \a 3 ), the same expressions (1331) - ([37D remain valid for 
the corresponding covariance of the operators a) (f), a 3 (t ): 


»)), - <a*(S)>„<%(«)>„ = (W))o = £ \V lt (t )\ 2 , (38) 

k 

where {0) a = («ia 2 |0|aia 2 ). Eq. (1331) is then independent of a.\ and a 2 . 

Similarly, we may evaluate the coordinates and momenta fluctuations and their 
dimensionless ratios to their initial values, 

M = «?(*)>« - WV)) 2 r2 ,, = (p 2 M)„ - (q(t)>; 

0)U «J(0))« - (<M0))r p ‘ u (if(o))„ - (P,(0)>J ’ 1 J 

which for a coherent initial state satisfy RQ j {t)Rp.(t) > 1, due to the uncertainty 
principle and the fact that a coherent initial state has minimum uncertainty. Squeezing 
in Qj or P 3 occurs whenever Rq.it) or Rp.(t)) becomes smaller than 1. We obtain, 
explicitly, 


R Qi(Pi){t) = 1 + 2[(if f (t)>. - l(0j(i))o| 2 ± Re«Oj(<))« - (aft)) 2 )] 

= 1 + 2[<A7;(t)) 0 ± Re((aJ(t)> 0 )], (40) 

where Re denotes real part, + (—) corresponds to Q 3 ( Pj ) and 

<<#)><> = £ • (41) 

k 

These ratios are then also independent of aq,a 2 , and deviate from 1 unless V{t) = 0. 

An initial coherent state is a pure separable gaussian state, which under the present 
Hamiltonian will remain gaussian (but no longer separable) Vf. Its entanglement entropy 
S(t ) can then be evaluated through the gaussian state formalism [23, 25] [38:, 39j[l0] and 
can be written as 


S(t) = -Tr P j(t) In p 3 (f) = - f(t ) In/(f) + [1 + /(f)] ln[l + /(f)] , (42) 

where pj{t ) denotes the reduced state of one of the modes and /(f) is the symplectic 
eigenvalue of the single mode covariance matrix: 

m = /««■,(())» - i(%(())«i 2 +1) 2 - !<<#)>.» - (“jM>ji 2 -1 

= v£(£££Fh( 7 w)T-i ( 43 ) 
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which is non-negative and the same for j — 1 or j — 2 if the state is pure and gaussian. 

It represents the effective occupation number of the mode [25]. It is obviously also 
independent of ai,a 2 , he., the same for any coherent initial state. The entanglement 
entropy (142|) is just an increasing concave function of fit). Again, in the isotropic case 
oq — uj 2 , V(t) — 0 (Eq. (J23])), entailing no generated entanglement when starting from 
I ol\ol 2 ). 


3.2. Results 



0 50 100 150 


Wit 



Figure 1. (Color online) Top: Temporal evolution of the entanglement entropy S(t) 
(l42l) and the shifted squeezing ratios (l44l) for the operators Qj, Pj of each oscillator, 
starting from a separable coherent state. We have set here u} 2 = uq/2 with a 
rotation frequency w = 0.49wi, such that the system is close to the first instability 
(occurring when u> reaches W 2 ). Both Q 1 and P 2 exhibit appreciable squeezing 
(A Oj(t) < 0), whose evolution is in phase with that of entanglement. Bottom: The 
corresponding symplectic eigenvalue /(f) (l43l) determining the entanglement entropy, 
and the average boson numbers of each mode when starting from the separable vacuum 
(or the covariances (l38l) when starting from a coherent state). Quantities plotted are 
dimensionless. 

Results for the previous quantities are depicted in Figs. [l]-[2] for u> 2 < uq. We 
concentrate on the Erst dynamically stable sector 00 < uj 2 . For improved visualization 
of squeezing, we use there the quantities 

A Qj(t) = R Qj (t) - 1, A Pj(t) = R P .(t ) - 1, 

with squeezing in Qj (Pj) indicated by a negative value of A Qj (A Pj). 


(44) 
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In Fig. |T]we consider the anisotropic case ca 2 = uq/2 with ca close to oa 2 , i.e., to the 
first instability, such that ca + ~ 1.31cui, ca- ~ 0.07u;i. The evolution then exhibits large 
amplitude low frequency oscillations governed by cu_, together with small amplitude high 
frequency oscillations governed by oa + . The picture clearly shows that in this regime 
squeezing and entanglement oscillate in phase : Maximum entanglement occurs at times 
t n ~ nnf(2ca-), n odd (see Eqs. (I45|l - (l46|l ). simultaneously with maximum squeezing 
in the operators Q\ and P 2 , and maximum average boson number in the oscillators (or 
maximum covariance (|38li in the case of an initially coherent state). 

These results can be approximately described by conserving just the main terms 
in Vjk(t) and Ujk(t ) (Eqs. (1T5h — (TT8li ) for small ca_, which are those proportional to cal 1 . 
We obtain 

w(o> o« , ( 45) 

<«#<)>o « - (i - . 

(46) 

where 8j is given in (TT9l) . Hence, both (Nj(t )) 0 and {a 2 (t ))o become proportional to 
sin 2 ca_t, entailing that all quantities plotted in Fig. [1] will be governed by this term, thus 
oscillating in phase. The presence of a factor caj 2 in (|45l) also entails (N 2 (t )) 0 > ( Ni(t )) 0 , 
i.e., the boson number will be larger in the mode with the smallest frequency, as verified 
in the bottom panel of Fig. [D 

Moreover, for co below but close to ca 2: the bracket in ()46l) will be positive for 
j = 1,2 and larger for j — 2 (since 0 < <5 2 < < 1), leading to (af(t)) 0 < 0 and 

(a 2 (t))o > 0, with |(a 2 (f)) 0 | < | (a| (t))o I - Hence, according to Eq. (HUD , squeezing will 
occur in Qi and P 2 , as seen in Fig. [H being more pronounced in P 2 . Eqs. (I45l) - (l46l) also 
show that squeezing and entanglement are driven by the anisotropy, i.e., they vanish for 
uj\ = co 2 and at fixed ca become larger as coi —ca 2 increases (or the ratio co 2 /loi decreases). 
They also increase, of course, as ca approaches ca 2 , i.e., as the first instability is reached 
(ca- = 0 for ca = ca 2 ). We also mention that for ca above but close to ca\ (i.e., in the 
second stable region but close to instability), the behavior is similar although squeezing 
will occur for P\ and Q- 2 , since the bracket in (|46l) becomes negative (with 5 2 > <5i > 1), 
implying (a?(f)) 0 >0, (a|(i)) 0 < 0. 

In contrast, away from instability (ca well below ca 2 ), quantum effects become much 
smaller even though they remain non-zero, as seen in Fig. [2] for ca = 0.15tai (where 
ua + ~ 1.04cji, ca- ~ 0.45o;i). Alternating squeezing in both Pj and Qj is now observed 
(the behavior of A P 2 and AQ 2 is analogous) and the correspondence with the evolution 
of entanglement (i.e. with f(t )) is less direct, with the maxima of f(t ) reflecting 
essentially the largest squeezing (that of Q i or P 2 ). Nonetheless, the average boson 
numbers ( Nj(t)) follow approximately f(t ). 

We can easily understand these results by considering the expansion of the exact 
expressions (|36l) - (|4ll) for small ca (|u;| Mirifwi, ca 2 ]). We obtain, neglecting terms of 
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t 


Figure 2. (Color online) Top: The shifted squeezing ratios for Q i and Pi for weak 
coupling ui = 0.15u;i. Now both Qj and Pj exhibit small alternating squeezing. 
Bottom: The corresponding value of f(t) (l43l) . together with the average boson number 
of each oscillator. The behavior of the entanglement entropy S(t) (l42l) is similar to 
that of /(f). 


order a; 4 , 


{NM c 


u 2 (ui - u 2 f 


sin 


2 UJ1+UJ2 


t, 


«?(*)> c 


UlL0 2 (uJi + UJ 2 ) 2 


— U)2 


uj-t l-UJo , 

sin — 1 
011+012 


(47) 
. (48) 


which show that both ( Nj(t)) 0 and (a|(t)) 0 are of order ta 2 , with (iV 1 (t)) 0 = (N 2 (t )) 0 at 
this order. This implies that |(a 2 (t))o| 2 will be of order ca 4 , so that Eq. leads to 
f(t) •ze (Nj(t)) 0 up to 0(u} 2 ). Hence, the ( Nj(t)) 0 and f(t) will be close for small u j. 

Besides, Re[(a 2 (t))o] will change its sign as t evolves, indicating that squeezing 
will alternate between Qj and Pj, being again larger for the oscillator with the lowest 
frequency due to the factor coj 1 in (l48]h Note that for small (Nj(t)) 0 and ( aj(t )) 0 , 
~ 1 + {Nj(t ))o ± Re[(aj(t) 2 )o], which will not strictly follow (Nj(t )) 0 , as 
(a 2 (f)) 0 is of the same order as ( N 2 (t)) 0 but not proportional to it. Finally, it is verified 
from ifTPi iflbl) that (Nj(t))o and ( a 2 (t)) 0 are again proportional to (cji — co 2 ) 2 and 
uj\ — respectively, hence vanishing for uji = oo 2 and leading to a larger entanglement 
and squeezing as the anisotropy uq — u 2 increases. 

We finally mention that for short times t such that oj±t 1, we obtain, after an 
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expansion np to 0(t 2 

<jy,(0>< 


of the exact expressions (I55D - (FHD . 
ui 2 (— U 2 )' 2 1 2 


4l0iUJ2 


(»m 0«i-i ) J+1 


OJ (OJ 


1 


w3) t 2 


4o;io;2 


(49) 


which are in agreement with the t — > 0 limits of Eqs. (I47ll - (l48l) . Hence, these quantities 
increase initially quadratically with time t, with (Ni(t))o = (N 2 (t))o ~ /(f) in this limit. 
Moreover, in this regime (af(t))o = — (a|(t))o is real, and positive if uj\ > 002 - This 
entails that squeezing will initially start in P\ (as seen in Fig. [2]) and Q 2 , with 


R Pl (t) « R Q2 {t) 


^ _ o; 2 (a;i — ^2)^2 ,2 
2u>\U)2 


(50) 


4. Conclusions 


We have derived the exact analytical closed form solution for the field operators of two 
linear oscillators coupled through angular momentum. We then applied the solution 
to investigate the relation between squeezing and entanglement generation in this 
model when starting from a separable coherent state. In the vicinity of instability, 
the generated entanglement between the modes shows a large amplitude-low frequency 
behavior (almost periodic), which is reflected in a similar behavior of the squeezing 
in the coordinate of one of modes and the momentum of the other mode. A different 
behavior occurs in the weak coupling regime, away from instability, where the generated 
entanglement is small and the squeezing is weak, exhibiting an essentially alternating 
behavior for the coordinate and momentum of each oscillator. Approximate analytical 
expressions describing these two regimes have also been derived from the general exact 
solution. 

The present solution has of course potential applications for studies of other 
quantum statistical properties such as higher order squeezing, antibunching of photons 
and other nonclassical photon statistics. The solution is also of interest for quantum 
information applications. As stated in the introduction, the present model admits 
distinct physical realizations, so that results could in principle be tested in quite different 
scenarios (optical simulations, particles in anisotropic harmonic traps, condensates, 
etc.). We remark, finally, that expressions similar to (TTOl) — (TTT1) and (13611 . (FXTjl remain 
formally valid for general systems of n harmonic modes interacting through quadratic 
(in aj, Oj) couplings, replacing % by the corresponding 2n x 2 n matrix. 
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Appendix 


By means of the canonical transformation (j = 1, 2) 

P'i = Pj + 7Q3-V Q'j = (Qj - v p 3-j )/(i + 'in 


(51) 


where 7 = 
m = 1 ) 


2A— -\-UJ2 


4 OJ 


, 77 = ■> we may rewrite ( 0 Q) in the decoupled form (we set here 


H = 


£ (« wf+Wj 2 ) . 


(52) 


3 = 1,2 


where ctj = 1 — ^(7 + (—lpu;), /3j = ^(7 — (—l) J 'o;). Here aq > 0, f3j > 0 for j — 1,2 

in the fully stable region a; < (n 2 , whereas aq > 0 , /?i > 0 , a 2 < 0 , fd 2 < 0 in the 

second dynamically stable sector u> > u>i, with CK 2 > 0, f3 2 < 0 in the unstable sector 

co 2 < to < oj\ and (3 2 = 0 at the borders uj — uj 2 or oj — uj\. Eq. (l52ll then leads, in the 

dynamically stable regions oj < co 2 or u> > uq , to 


Pi + Ql Pl + Ql 
H = huj+ + —— ± hjj_- 


(53) 


2 2 

where oj± = (real), P± = a 1 / h(3i Pi , Q± = ^J]h/iionQi, and the minus 

sign in (|53|) applies for to > oj\. Eqs. (T33l) (T34|) are then obviously obtained for 
b± = (Q± + iP±)/\/2, b± = (Q± — iP±)/V 2 . The possible normal representations 
in the unstable regime are discussed in detail in 
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